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ABSTRACT 



Context. Cygnus X-3 (Cyg X-3) is a well-known microquasar producing variable emission at all wavelengths. Cyg X-3 is a prominent 
X-ray binary producing relativistic jets, and studying its high energy emission is crucial for the understanding of the fundamental 
acceleration processes in accreting compact objects. 

Aims. Our goal is to study extreme particle acceleration and y-ray production above 100 MeV during special spectral states of Cyg X- 
3 usually characterized by a low hard X-ray flux and enhanced soft X-ray states. 

Methods. We observed Cyg X-3 with the AGILE satellite in extended time intervals from 2009 Jun.-Jul, and 2009 Nov.-2010 Jul. We 
report here the results of the AGILE y-ray monitoring of Cyg X-3 as well as the results from extensive multiwavelength campaigns 
involving radio (RATAN-600, AMI-LA and Metsahovi Radio Observatories) and X-ray monitoring data (XTE and Swift). 
Results. We detect a series of repeated y-ray flaring activity from Cyg X-3 that correlate with the soft X-ray states and episodes of 
decreasing or non-detectable hard X-ray emission. Furthermore, we detect y-ray enhanced emission that tends to be associated with 
radio flares greater than 1 Jy at 15 GHz, confirming a trend already detected in previous observations. The source remained active 
above 100 MeV for an extended period of time (almost 1.5 months in 2009 Jun.-Jul and 1 month in 2010 May). We study in detail 
the short timescale y-ray flares that occurred before or near the radio peaks. 

Conclusions. Our results confirm the transient nature of the extreme particle acceleration from the microquasar Cyg X-3. A series of 
repeated y-ray flares shows correlations with radio and X-ray emission confirming a well established trend of emission. We compare 
our results with Fermi-LAT and MAGIC TeV observations of Cyg X-3. 
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1. Introduction 

Cygnus X-3 (Cyg X-3) is one of the most puzzling and 
interesting com pact objects i n our Galaxy. Discovered 
in 1966 (.Giacconi et al. 1967h . Cyg X-3 is a high- 
mass X-ray binary exhibiting a 4.8 hr modulati on in 
its X-ray (iParsign ault et al. 1972t IVilh u et al 2009h. in - 
frared _^ason et al. 19861; Ivan K erkwii k et al. 1992 



Ivan Kerkwiik et al. 199*61) and y-ray (lAbdo et al. 2009 
emission. The system is believed to be composed of a 
mass-donating Wolf-Rayet star (Ivan Kerkwiik et al. 19961) 
orbiting around a compact object (neutron star or black 
hole). Certainly what makes Cyg X-3 an interesting ob- 
ject is its abili ty to produce v ery energetic relativistic 
radio jets (Molnar et al. 19881: iMiodu szewski et al. 2001) 
that us ually occur in coincid e nce with particular spectral 
states jMcColl ough et al. 19991: IZdziarski & Gierlinski 2004*: 



Szostek et al. 2009 



Koljonen et al. 2010). Many 



gghetal. 1999t |z^ 
'10). Many decade 



Hialmarsdotter et al. 2009: 
decades of X-ray monitoring 
of Cyg X-3 have provided a wealth of information on this 
source, and in general the soft, intermediate and hard X-ray 
states show correlated or anti-correlated behaviours with respect 
to the radio and jet emission. In particular, the soft (1-10 keV) 
and hard (20-lOOkeV) X-ray emission from Cyg X-3 are clearly 
anti-correlated during "normal" stages during which the source 
produces radio emission at low or intermediate levels with no 
major radio jet production. Major radio flares in Cyg X-3 are 
preceded by quenched states, during which the source is in a 
soft X-ray state and radio emission is strongly suppressed. 

The ability of Cyg X-3 to efficiently accelerate particles 
in relativistic jets, and the favourable jet geometry parame- 
ters that make the transient jet phenomenon very dramatic (e.g. 
(Mioduszewski et al. 2001) ) have attracted considerable atten- 
tion from the astrophysics community. Cyg X-3 is a Galactic 
"microquasar", and the inner dynamics and accretion properties 
of its compact object can shed light on fundamental physical pro- 
cesses of relativistic jet sources, both Galactic and extragalac- 
tic. From this point of view, exploring the high-energy emission 
from Cyg X-3 and unveiling its temporal and physical properties 
open the way to a detailed understanding of the plasma proper- 
ties of inner accretion disks of compact objects and of particle 
acceleration processes under extreme radiative conditions. 

Motivated by these reasons, our group has made a detailed 
study of y-ray emission from Cyg X-3 with the AGILE satellite. 
The AGILE discovery of transient y-ray emission from Cyg X-3 
in 2008 Apr associated with a specific spectral state preceding 
a major radio jet ejection opened a new window of investigation 
of microquasars. Several other major y-ray emission episodes 
from Cyg X-3 have been detected by AGIL E and Fermi since 

2008 (i Tavani et al. 2009bl: lAbdo et al. 2009b'). In this paper, we 
focus on the main results of our extensive search for transient y- 
ray emission from Cyg X-3 carried out in the energy range 100 
MeV - 50 GeV by AGILE during the periods 2009 Jun-Jul and 

2009 Dec.-2010 mid-Jun. We find that the activity in y-rays dur- 
ing the periods 2009 Jun-Jul and 2010 May temporally coincide 
with the hard X-ray emission reach minimum values. O ur results 
confirm the overall trend found by (iTavani et al. 2009bh . i.e., that 
the highest y-ray emission above 100 MeV from Cyg X-3 is as- 
sociated with soft X-ray spectral states that are in general co- 
incident with or anticipate radio jet ejections. Furthermore, our 
results show that y-ray emission from Cyg X-3 is detectable by 
AGILE not only during relatively short (1-2 day) flares as in 



(ITavani et al. 2009bl) but also during "extended" periods lasting 
several days or weeks (as during 2009 Jun-Jul). Detecting con- 
tinuous y-ray emission during "active" phases is of great theo- 
retical relevance for the modelling of Cyg X-3 . 

We also briefly compare here the AGILE and Ferm/-LAT 
results. We also consider a recent investigation by the MAGIC 
group reporting their results on a monitoring program of Cyg X- 
3 in the TeV energy range (Aleksi c et al. 2010), and briefly dis- 
cuss the implications with respect to the AGILE detections. 

2. The AGILE GRID observations of Cyg X-3 

AGILE (Astrorivelatore Gamma ad Immagini LEggero - Light 
Imager for Gamma-ray Astrophysics) is a scientific mission of 
the Itahan Sp ace Agency (ASI) launched on April 23, 2007 
(ITavani et al. 2009a). The AGILE scientific payload is made 
of three detectors: (1 ) a y-ray imager made of a Tungsten- 
Silic on Tracker (ST) dBarbiellini et al. 20021: iPrest etal. 200l 
Bul garelh et al. 201 Oh with a large field of view (about 60°); 
(2) a co-axi al hard X -ray silicon detector (named Super- 
AGILE (Fer oci et al. 20 07)). for imaging in the 18-60 keV en- 
ergy ran ge, and (3) a C sI(Tl) Mini-Calorimeter (MCAL) de- 
tector CLabanti et al. 2006i) that detects y-rays or particle en- 
ergy deposits between ~350 keV and 100 MeV. The whole 
instrument is surrounded by an anti-coincidence (AC) system 
(Perotti et al. 2006) of plastic scintillators for the rejection of 
background charged particles. An effective background rejec- 
tion, event trigger logic, and on-boar d data storage and trans- 
mission is implemented (lArgan et al. 2 004). ST, MCAL and AC 
form the so called Gamma-Ray Imaging Detector (GRID) for 
observations in the 30 Me V-50 GeV y energy range. The AGILE 
orbital characteristics (quasi-equatorial with an inclination an- 
gle of 2.5 degrees and average 530 km altitude) are optimal for 
low -background y-ray observations. AGILE data are transmitted 
to the ASI Malindi ground station in Kenya, and quickly trans- 
ferred to the ASI Science Data Center (ASDC) near Frascati. 
Data processing of y-ray data is then carried out at the ASDC 
and AGILE Team locations. 

The AGILE-GRID is optimized in the 100 MeV - 
1 GeV range, as demonstrate d by the calibrated data 
(ICattaneo et al. 20Tl1; Chen et al. 201 lb) and by the in-orbit per- 
formance (ITavani et al. 2009al) . 

The AGILE y-ray exposure of the Cygnus region can be di- 
vided in two parts. 

(1) A first set of exposures obtained in the satellite "point- 
ing mode" characterized by pointing centroids near the Cygnus 
region in the Galactic plane. In this mode, AGILE accumulated 
a total effective time of ~338 days during the period 2007 Jul- 
2009 Oct. Table [llprovides the details of the main time intervals 
analyzed in this paper for which a substantial exposure was ac- 
cumulated in the Cygnus region in this pointing mode. 

(2) A second set of exposures obtained quasi-continuously 
with the satellite operating in "spinning mode" since 2009 Nov. 
In this mode, the satellite axis sweeps an entire circle in the sky 
in approximately 7 minutes. Depending on the season, the whole 
sky is progressively exposed with a typical accumulating pat- 
tern. The Cygnus region is favorably positioned in the sky, and 
has been continuously monitored since 2009 Nov. with a 2-days 
exposure comparable with 1-day pointing exposure level. 

2.1. Data analysis method 



Send offprint requests to: A. Bulgarelli, e-mail: 
bulgarelli@iasfbo . inaf . it 



The data analysis was performed on the data set in pointing 
mode generated with the reprocessing no. 3 of the AGILE 
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Table 1. AGILE observations of the Cyg X-3 region analyzed in this paper. 



Obs. block 




b 


Time (MJD) 


Mean off-axis angle 




Flux 


OB7500/OB7600 


92.835 


-9.574 


54997.50-55027.50 


17/29 


3.5 


32 ± 10 


OB7700 


105.729 


7.231 


55027.50-55055.50 


30 


<2 


<20 


spinning 






55168.50-55362.50 




<2 


< 16 



The Table provides: (1) the observation block (OB) number; (2-3) the galactic coordinates / and b of the pointing centroids; (4) the time interval 
in MJD; (5) the off-axis position of Cyg X-3 at the beginning of the OB (in degrees); (6) the statistical significance of the source detection 
according to the maximum likelihood ratio test; (7) the period-averag ed flux F {E > 100 MeV) in lO"** ph. cm^^ s-'(if ^JT, < 3 a 2 - o" upper 
limit is reported) 



Standard Pipeline, and with the AGILE-GRID software pack- 
age version 4 publicly available at the ASI Data Center web site 
(http://agile.asdc.asi.it/). The analysis has been performed with 
the FM3.119 filter and the calibration matrix used is the 10023 
version. The events collected during the passage in the South- 
Atlantic Anomaly and the Earth albedo background were con- 
sistently rejected. The GRID event direc tions were reconstru cted 
by the AGILE Kalman filter technique jGiuliani et al. 20061) . To 
reduce the particle background contamination, we selected only 
events flagged as confirmed y-ray events (G class events, corre- 
sponding to an on-axis effective area of ~ 350 cm^ at 100 MeV). 

The multi-sour ce likelihood analysis method (MSLA) 
dMattox et al. 1996h was used to search for persistent and tran- 
sient emission from Cygnus region; this analysis method iter- 
atively optimizes position, flux and spectral index of all the 
sources of the region. The likelihood ratio is then simply 
the ratio of the maximum likelihood of two hypothesis (e.g. the 
absence and the presence of a source). For this analysis, inte- 
grating all flaring episodes in 2009 Jun-Jul, we obtained (see 
Section [22] and Fig.|4| and kept fixed the photon index of 2.0 
for Cyg X-3 . This is consi stent with other AGILE detections 
reported in dTavani et al. 2009b) . During the analysis the posi- 
tion of the source is kept free and constrained with the 95% 
confidence contour level. The AGILE photon counts, exposure, 
and Galactic background maps were generated with a bin size 
of 0.3° X 0.3° for E>100 MeV to compute the period-averaged 
source flux and its evolution. The analysis was performed over a 
region of 10° radius. 

The Galactic diffuse y-ray radiation jGiuliani et al. 20d4h 
and the isotropic emission are taken into account in the model. 
The Galactic diffuse emission model is based on a 3D grid with 
bin of 0.25° in galactic longitude and latitude and 0.2 kpc in 
distance along the line of sight. In order to mode l the matter dis- 
tribu tion in the galaxy we use the HI (ikalberla et al. 20051) and 
CO jDame et al. 200 lb radio survey an d we take into account 
the cosmic ray models (IChi et al. 199lh which can differ from 
the locally observed cosmic ray spectrum. 

Two energy bands have been considered in the analysis: 100 
MeV - 50 GeV and 400 MeV - 50 GeV 

Only the detections with -yTJ ^ 3 and with a position con- 
sistent with Cyg X-3 source are selected. For this subclass of 
selected flare peaks, a non-automatic verification is performed 
to further confirm the results. Pre-trial significance was deter- 
mined b y calculating the density function of our analysis 
method (Mattoxet al. 1996) by means of Monte Carlo simula- 
tion (Bulgarel li et al. submitted) . For the determination of the 
post-trial significance we have taken into account only the num- 
ber of bins of the light curv43- 



' This means that we have considered 58 maps in pointing mode and 
92 maps in spinning mode. The total number of maps is 150. 



The Cygnus region is characterized by complex y-ray emis- 
sion: several y-ray sources are detected above 100 MeV, and it is 
important to correctly model the diffuse radiation of the region. 
Three bright y-ray sources dominate the Cygnus region and also 
reported in the First A GILE Cata log (Pittori et al. 2009). They 
are all y-ray pulsars ( iHalpern et al. 2008i ; lAbdo et al. 2009at 
ICamilo et al. 2009). We added more nearby sources in the y-ray 
model (Table |2]l and the statistical analysis of these sources is 
compared with the First AGILE Catalog analysis. A new likeli- 
hood analysis investigation has been performed to take into ac- 
count a new background event filter (FM3.119). We have also 
added considerably more exposure to the data set. 

For our MSLA method we first determined self-consistently 
the positions of all bright sources (including also a Cyg X-3 can- 
didate source, AGL 2033+4056) by considering the whole set 
of AGILE-GRID data in pointing mode for the Cygnus region. 
Subsequently (see next Section), for the time-resolved analysis 
of Cyg X-3 candidate source, we fixed the fluxes and the posi- 
tions of the brightest sources as shown in Table |2] keeping the 
position and flux of Cyg X-3 free. In particular, we studied the y- 
ray sourceQthat we call here AGL 2030-H4129 that is positioned 
only 0.4° from Cyg X-3 position. This y-ray source is consistent 
in position and ave rage flux with the source IFGL J2032.2-H4127 
('Ab do etaL 20 10*) which is the counterpart of the radio pulsar 
PSR J2032-I-4127 (of sp in period P = 147 ms, and a most likely 
distance d ^ 1.8 kpc dCamilo et al. 2 009)). Given the differ- 
ences in the published determination^ of the y-ray spectral in- 
dex of AGL 2030+4129/lFGL J2032.2-H4127, we considered in 
our analysis two y-ray spectral indices - 2.24 (Ferm/-LAT 
First Catalog) , and a2 - 1.84 + 0.2 calculated with AGILE data; 
for these indices we find the corresponding average fluxes 26 + 3 
and 18 ±4 ■ 10"^ ph. cm"^ s"' . In this paper we adopt the spectral 
index value a = 1.84 + 0.2 that is the value calculated with the 
AGILE data. 

An alternative method using the False Discovery Rate (FDR , 
(iBeniamini et al. 19951: iMiller et al. 20"oTI: iHopkins et al. 2002h ) 
has been used to analyze the flares detected with MSLA. The 
detection method (FDRM) is a statistical test that takes into ac- 
count the corrections for multiple testing, as needed for example 
in repeated systematic searches. The FDRM allows to control 



' Notice that the initially re ported AGILE First Catalog source 
lAGL J2032+4102 jPittori et al. 2009) comprises two y-r ay sources 
separated by only 0.4° as clarified in dTavani et al. 2009bh and con- 
firmed by our more refined analysis. 

Fermi-hAT First Catalog analyses give a y-ray photon spectral in- 
dex consistent with the value a = 2.24 ± 0.034. However, an analysis of 
the y-ray Fermi-LAT data of PSR J2032+4127 that was carried out by 
selecting only the on-pulse phase photons gives an exponentially cut- 
off power law spectrum of photon spectral index a = 1.1 ± 0.2 ± 0.2 
and cutoff energy E^ = (3.0 ± 0.6 ± 0.7) GeV for a pulsed y-ray flux o f 
= (7 ± 1 ± 1) X 10^'* cm-' S-' above 100 MeV dCamilo et al. 20091) . 
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Table 2. List of Cygnus region sources with E > 100 MeV. 









E>100 MeV 






E>400 MeV 








AGlLt Name 


1 


b 




Flux 






Flux 


Counterpart name 




AGL 2021+4029 


78.24 


2.16 


42.1 


141 : 


t 4 


30 


33 ± 1.5 


lAGL J2022+4032 - 


Gamma Cygni 


AGL 2021+3652 


75.24 


0.14 


23.3 


67 ± 


3 


21 


19 ± 1.2 


lAGL J2021+3652 - 


PSR J202 1+3651 


AGL 2030+4129 


80.11 


1.25 


8.1 


18 ± 


4 


8.7 


5.6 ± 0.7 


lAGL J2032+4102 - 


LAT PSR J2032+4127 


AGL 2026+3346 


73.28 


-2.49 


6.8 


10 ± 


1.7 


6.2 


2.8 ± 0.5 






AGL 2046+5032 


88.99 


4.54 


6.5 


10 ± 


1.7 


6.1 


2.6 ± 0.5 












O.J 


1 A -4- 


Z.J 










AGL 2029+4403 


81.97 


3.04 


5.4 


14 ± 


3 


5.1 


3.4 ± 0.7 






AGL 2033+4056 


79.92 


0.58 


5.3 


15 ± 


2 


3.1 


2.3 ± 0.7 


lAGL J2032+4102 - 


Cygnus X-3 


AGL 2038+4313 


82.32 


1.18 


5.1 


15 ± 


3 










AGL 2024+4027 


78.56 


1.63 


5.0 


24 ± 


5 


6.1 


11 ±2 






AGL 2019+3816 


76.24 


1.14 


4.2 


11 ± 


2.4 


4.2 


2 ±0.6 






AGL 2036+3954 


79.47 


-0.56 


3.4 


5.0 4 


: 1.5 


5.7 


3.8 ± 0.7 







The Table provides: (1) AGILE name of the sources; (2) (3) the galactic coordinates / and h; (4) (6) the statistical significance s/T^ of the source 
detection according to the maximum likelihood ratio test for E>100 MeV and E>400 MeV respectively; (5) (7) the period-averaged flux F 
(E > 100 MeV) in 10"** ph. cm"- s^'for E>100 MeV and E>400 MeV respectively; (8) a possible counterpart 



the expected rate of false detections (due to background fluc- 
tuations) within a selected sample. The FDR-a parameter pro- 
vides the fraction of expected false detections for a given source 
selection. The method was adapted to the analysi s of AGILE 
gamma-ray data of the Galactic plane (Sabatini et al. 20l"ol: 
iTavani et al. 2 009b) in pointing mode and for E>100 MeV (for 
spinning mode end for E>400 MeV the method is not avail- 
able). The FDRM is complementary and more conservative than 
the determination of a post-trial detection significance based 
on a simulated set of replicated photon maps; applying FDR 
method to typical short-time AGILE-GRID photon maps in the 
Galactic plane, FDR-o' values near or below 0.01 fully qual- 
ify for gamma-ray transients (it correspond to post-trial random 
occuiTences of 1-day map replications equal to 1 every 300 or 
more). 



2.2. Flaring y-ray activity observed by AGILE 

We detected several significant short timescale y-ray flares from 
the Cyg X-3 region during the 2009 Jun-Jul period (pointing 
mode) and during the 2009 Dec. -2010 Jun period (spinning 
mode). Fig. [1] shows the light-curve of Cyg X-3 obtained in the 
15-50 keV range by the Swift -BAT instrument from 2008 Jan. 1 
to 2010 Jun 1. Superimposed on the same plot are several arrows 
indicating the periods of major y-ray flaring detected by AGILE 
an d associated with Cy g X-3, including those already reported 
in (ITavani et al. 2009bl) (see Table|5ll. 

Table|6]and|7]reports the details of the flares analyzed in this 
paper with the significance of the detections both with MSLA 
method and FDR method (when applicable). The MSLA and the 
FDRM use different assumptions for the background (diffuse) 
model, therefore it could happen that the relative significance 
of the detections is different. However, the result we wanted to 
point out in the table is that both independent methods would 
have selected the episodes as statistically significant. 

Periods of long y-ray AGILE exposure (shown in grey in 
Fig. [TJ cover both high, intermediate, and low states of hard X- 
ray emission from Cyg X-3. 



2.3. Average y-ray flux detected from Cyg X-3 

By integrating all AGILE data from 2007 Jul 13 (MJD 
54294.50)to 2010 Jun 15 (MJD 55362.50) (about 1.5 year of ef- 
fective Uvetime on the Cygnus region) for E> 100 MeV. we find 
a y-ray source at the ^/^\ - 5.3 level consistent with Cyg X- 
3 position, at (l,b) = (79.92, 0.58) + 0.3 (stat.) ± 0.1 (syst.) and 
with an average flux F[ = (15 ± 2) ■ 10"'^ ph. cm"^ s^'above 100 
MeV that contains all the y-ray states. No significant detection 
is made for E> 400 MeV. 

If we remove all flares d etected by AGILE with V^T > 3 
(reported in ( Tavani et al. 2009b) and in Table |6l) we find a y-ray 
source at the = 4.5 level at ±e position (l,b) = (80.03, 0.65) 
± 0.5 (stat.) + 0.1 (syst.), consistent with the Cyg X-3 position, 
with an average fluxF2 = (11.3 ± 2.9) ■ 10"^ ph. cm"^ s^'above 
100 MeV. Note that this extended time interval includes sev- 
eral periods with low hard X-ray/high soft X-ray emission that 
we call here "active y-ray states". It is therefore possible that 
the value of F2 is influenced by low-level y-ray emission from 
Cyg X-3 during the active states. 

Table |5]reports all the "active" and "non-active" y-ray states 
associated with the hard and soft X-rays state. These periods 
are also shown in the bottom panel of Figure [1] For example, 
integrating during the period 2009 Jun 15 until 2009 July 15 
(MJD = 54997.50 - 55027.50) we find an average flux F3 = 
(32 + 10) ■ lO-*^ ph. cm-2 s"' above 100 MeV. 

In order to study "non-active y-ray" Cyg X-3 states (i.e., 
those corresponding in general to the common high hard X- 
ray states) we also focused on specific time intervals charac- 
terized by average Swift/BAT levels of emission above 0.028 
cts. cm"^ s For the period starting on 2009 Jan. 19 until 2009 
Feb. 28 (MJD = 54850.50-54890.50) we find a 95% confi- 
dence limit upper limit F4 < 23 ■ 10"^ ph. cm"^ s"' above 100 
MeV. Integrating from 2009 Dec. 1 until 2010 May 4 (MJD = 
55166.50-55320.50) we find a 95% confidence limit upper limit 
F5 < 15 • 10-** ph. cm-2 s-i above 100 MeV. 

2.4. Instrument stability 

We have analyzed some AGILE sources having flux near the av- 
erage y-ray flux detected from Cyg X-3 (see Section l231 l that are 
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Time [MJD] 

Fig. 1. (Top panel) Swift /BAT lightcurve ( in counts per second in the energy range 15-50 keV) and AGILE-GRID y-ray flares for E> 
100 MeV as a function of time. The red arrows mark the dates of major y-ray flares ofCyg X-3 : MJD 54572, 54772, 54811, 54998, 
55001, 55003, 55007, 55025, 55034, 55324, 55343. Gray areas show the interval of good AGILE y-ray exposure of Cyg X-3 . 
(Bottom panel) Average y-ray flux from Cyg X-3 region in hard and soft X-ray state for E>100 MeV. 



Table 3. Average y-ray flux (E> 100 MeV) from Cyg X-3 region 
in hard and soft X-ray state. 



Time (MJD) 


V7T 


Flux 




X-ray state 


54566.50-54646.50 


5.9 


36 ± 


7 


soft 


54770.50-54850.50 


3.6 


22 ± 


6 


soft 


54850.50-54890.50 


<2 


<23 




hard 


54936.50-54980.50 


<2 


< 19 




hard 


54997.50-55027.50 


3.5 


32 ± 


10 


soft 


55027.50-55055.50 


<2 


< 20 




hard 


55166.50-55320.50 


<2 


< 15 




hard 


55320.50-55350.50 


<2 


<59 




soft intermediate 


55350.50-55400.50 


<2 


<27 




hard 



The Table provides: (1) the time interval in MJD; (2) the statistical sig- 
nificance yfT\ of the source detection according to the maximum like- 
lihood ratio test; (3) the period-averaged flux F (E > 100 MeV) in 10"** 
ph. cm"^ s"'(if yfT\ < 3 a 2 - cr upper limit is reported); (4) the X-ray 
state 



not expected to be variable. In particular, we have considered the 
y-ray pulsars reported in Table H] The analysis has been carried 
out using the same procedure performed on Cyg X-3, by divid- 
ing the analyzed period (from 2007 Jul to 201 1 May) in two sets: 
pointing mode (with a 1 day bin size light-curve) and spinning 
mode (from 2009 Nov., with a 2 days bin size light-curve). We 
have done about 1530 trials for the pointing mode and 1536 trials 
for the spinning mode. The resulting probability density function 
(PDF) of values is shown in Fig. |2]where it is compared with 
a simulation of Cyg X-3 field: in this field all the sources re- 
ported in Table |2] are simulated except Cyg X-3 that has zero 
flux. It is noted that the distribution of real AGILE observa- 



tions of the source listed in Table|4](red line) is fully compatible 
with the simulations (black line) and this fact excludes the pres- 
ence of spurious flares introduced by in strument instability dur- 
ing the AGILE lifetime. As stated in (Bulgarelli et al . submittedh 
the vertical translation of the T, distribution depends on the com- 
plexity of the analyzed region. 

With these distributions the expected number of wrong de- 
tections with yfiTs) > 3.1 (the lowest value of Table|6]l in 150 
maps (the number of bins of the light curves) is about 0.2; in this 
paper we have 8 detections in 150 maps. 

2.5. Observations in pointing mode 

The AGILE pointing observations of the Cygnus field analyzed 
in this paper are listed in Table [1] and cover about two months 
of uninterrupted observation. Having obtained the general re- 
sults on the Cygnus region described in the previous section, 
we carried out a systematic search for short timescale variabil- 
ity of Cyg X-3 during the interval of substantial AGILE-GRID 
exposure in pointing mode in 2009. 

Fig. [3] panels 1-3 provide a close-up of Fig. [T] showing the 
details of the y-ray emission above 100 MeV (panel 1) and 
400 MeV (panel 2) detected by AGILE from Cyg X-3 together 
with simultaneous hard X-ray information from both Swift -BAT 
and Super-AGILE (panel 3). The interval covering the period 
2009 Jun.-Aug. (MJD: 54997.50 - 55055.50) corresponds to a 
deep and prolonged minimum of the hard X-ray emission from 
Cy g X-3. The Fermi Ught curve is also superimposed in panel 
1 (lAbdo et al. 2009bl) . Several 1-day episodes of y-ray flaring 
emission from Cyg X-3 with flux F > 100 x 10 ph. cm"^ s"' 
are detected. Table|6]provides detailed information on the main 
flare episodes of Fig.[3]and obtained for 1-day time integrations 
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Fig. 2. The red histogram is the probability density function 
(PDF) ofTs values for real AGILE observations of the sources 
listed in Tabled The black histogram is the PDF of the T^ val- 
ues of a simulated Cygnus region with all the sources reported in 
Table\2\except Cyg X-3 that has zero flux. The notch in the distri- 
bution near Ts—6 is caused by the switch between the fixed and 
the free position of the source in our analysis method. The red 
dotted line is the theoretical distribution, the green dashed 
line is the x\ theoretical distribution, the cyan dotted-dashed line 
is the jXj, distribution. 

for £ > 100 MeV. Table|7]lists flares for £ > 400 MeV. The c ov- 
erage of one of the MAGIC observations (l Aleksic et al. 2010h of 
Cyg X-3 in the soft X-ray state that provides an upper limit above 
250 GeV is also shown in panel 2. 

Fig. [3] panel 5-6 shows the radio monitoring data obtained 
by our group at 15 GHz at the AMI- LA radio telescope (al- 
ready published in (lAbdo et al. 2009bh ). the RATAN-600 data 
at the frequencies 2.15, 4.8, 7.7, 11.2, and 21.7 GHz, and 
the Metsahovi Radio Observatory data at 37 GHz. Finally, 
Fig. [3] panels 4-5 shows the 1.3-12 keV data obtained from the 
XTE/ASM monitoring of Cyg X-3. 

Fig. m shows the y-ray spectrum obtained by integrating all 
of the y-ray flaring episodes in 2009 Jun.-Jul. A single power- 
law fit gives a photon spectral index a - 2.0 + 0.2. However, a 
more complex spectrum with substantial curvature in the energy 
range 100 MeV - a few GeV cannot be excluded. 

2.6. Observations in "spinning" mode 

Since early 2009 Nov. the AGILE satellite has been operating 
in a "spinning" mode implying a smooth and continuous change 
of the satellite attitude. In this mode, the y-ray boresight axis 
sweeps 360 degrees in about seven minutes. Solar panels are 
kept perpendicular to the Sun direction by an automatic mech- 
anism, so that the pattern swept on the sky slowly moves with 
time following the solar panel configuration. The resulting y- 
ray daily exposure covers about 70% of the whole sky every 
day (leaving uncovered only the regions near the Sun or anti- 
Sun directions) and provides significant continuous monitoring 
of exposed sources for many months. The Cygnus region was in 
a favorable position in the sky, i.e., always in the field of view 
(within about 50 degrees off-axis) of the AGILE spinning in- 
strument. This was useful for assessing in a consistent way the 



pattern of y-ray variable emission from Cyg X-3 and correlating 
with other wavelength emission. The AGILE pointing mode had 
a relatively larger daily exposure compared with that obtained in 
the spinning mode. However, the AGILE pointings at the Cygnus 
region required a pre-defined planning strategy, and the overall 
monitoring livetime was about 50% from 2007 Jul to 2009 Oct. 

Since AGILE started operating in spinning mode (2009 
Nov.) to 2010 Jun. Cyg X-3 has been mostly in its hard 
X-ray state that does not favour strong y-ray emission 
dTavani et al. 2009bl:lAbdo et al. 2009bl) . 

We continuously searched for transient y-ray emission from 
the Cyg X-3 region during the AGILE spinning mode phase with 
an automatic alert monitoring system (BulgarelU et al. 2009i). 
In 2010 May this monitoring system detected a signal above 
Vfr.!) = 3.5 from Cyg X-3, that was subsequently ver- 
ified by a manual analysis. We detecte d two y-ray emis- 
sion episodes during MJD 55324-55326 (Bukar elh et al. 2010al) 
not confirmed by Fer mi (Corbel et al. 2010ai) , and on MJD 
55343-55345 (Bukar elli et al. 2010bl) confirmed by Fermi 
(ICorbel e t al. 2010b). The details of a refined analysis are re- 
ported in Table |6] for £ > 100 MeV and in Table [T] for E>400 
MeV. 

Fig. |5] shows the detailed 1-day hard X-ray lightcurve of 
Cyg X-3 together with the AGILE-GRID emission episodes and 
AMI-LA radio monitoring. A 15 GHz radio flare is evident in the 
data and it anticipates the y-ray flare at MJD 55343 by ~l-day. 

2. 7. Combining all Cygnus X-3 flares detected by AGILE 

Fig.|6]reports the sum of all the flares detected by AGILE-GRID 
from Cyg X-3 reported in this paper for E>100 MeV. The sig- 
nificance of the sum of all the flares is ^/^\ - 6.2 with a flux 
of 160 ± 40 ph. cm"^ s^'and a 95% contour level centered in 
(I, b) - (79.69, 0.72) in Galactic coordinates, with a semi-major 
axis of 0.47° and a semi-minor axis of 0.28° (including statistical 
and systematic errors). 

The significance of the sum of all the flares for E>400 MeV 
is yjTl - 5.1 with a flux of 58 ± 18 ph. cm"^ s^'and a 95% con- 
tour level centered in (I, b) - (79.8, 0.5) in Galactic coordinates, 
with a semi-major axis of 0.53° and a semi-minor axis of 0.34° 
(including statistical and systematic errors). 

The IFGL J2032.2+4127 and AGL 2030-1-4129 are both out- 
side the AGILE error box. This fact excludes that the detected 
flares originate from these nearby sources. 

In order to assess the statistical significance of our detec- 
tions, we consider the post-trial probability of flare occurrence. 
We have to distinguish two cases: 

1 . the case of a single flare episode originating from a specific 
source within a given error box (that we define as "single 
independent occurrence"); 

2. the case of repeated flaring episodes originating from a spe- 
cific source with a given error box (that we call here "re- 
peated flare occurrence". 

For each individual detection by AGILE reported in Tables|6]and 
|7]we calculated the post-trial significance of the single indepen- 
dent occurrences, which does not take into account the history 
of repeated occurrences. 

We calculate the post-trial significance for repeated flare oc- 
currences at the Cyg X-3 error-box position as follows. Each 
independent time period is a single trial. The chance probabil- 
ity of having k or more detections at a specific site with a T^ 
statistic satisfying > h in N trials is given by P{N,k) = 
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1 - Y/jJ) y j j P^^^ ^ where p is the p-value corresponding 

to the h value. For Ts > 10.9, we have a p-value of 6.8 x 10""*. 
The estimated probability of 8 detections (E>100 MeV) consis- 
tent with the null hypothesis in 150 maps is P = 2.3792 x lO"^-* 
that corresponds to ~ 7.2 Gaussian standard deviations. 

Table |2]reports the repeated flare occurrence significance for 
E>400 MeV. 

In order to test the pre-trial type I error (rejecting the null 
hypothesis when in fact it is true) we also searched for transient 
emission at the position (l,b)=(78.38, 0) near Cyg X-3 with sim- 
ilar characteristics of the diffuse background. We found a dis- 
tribution compatible with the expected statistics (the black his- 
togram of Fig. |2]i. 



Correlating hard X-rays and y-rays we find a negative cor- 
relation with a 3cr confidence level between these frequencies 
compatible with + 8 days lag. 

2.8.2. Correlation with hard X-ray emission 

Despite some gaps in the radio data, in particular during the 
pointing mode observations, several radio flares are evident in 
the data. We find a positive correlation above the 3cr confidence 
level between the radio and y-rays with a lag of + 8 days. The 
same consideration reported for the DCF between X-rays and 
y-rays is valid for the correlation between the radio and y-rays, 
and for this reason the DCF is not used for a measurement of the 
lag. 



2.8. Statistical correlations with other wavelenghts 

2.8.1. Anticorrelation with hard X-ray emission 

We notice that the y-ray episodes detected by AGILE are all in 
correspondence with peculiar states of the hard X-ray emission. 
Both the y-ray emission flaring episodes correspond to ~ 1-day 
minima of the BAT hard X-ray light-curve. The typical trend 
of anticorrelated soft and hard X-ray emission from Cyg X-3 is 
evident and also shown by the ASM hardness ratio trend. 

This paper reports eight flares originating from Cyg X-3 for 
E>100 MeV. All are in correspondence with a low hard X-ray 
flux of Cyg X-3; seven of them are under 0.02 cts. cm"^ s 'in 
the Swift/BAT, and one is under 0.028 cts. cm"^ s"' . To quantify 
the relation between y-ray emission and hard X-ray emission 
it is reasonable to set a threshold of 0.028 cts. cm"^ s^'in the 
Swift/BAT data: the fraction of the time that this source is under 
this level during the AGILE observation time is / = 0.32. For k 
y-ray flares originating from Cyg X-3, is the probability that 
all the flares are under the hard X-ray threshold, and P(N, k) is 
the probability of having k or more detections in maps (see 
Section lZTt . We find that the confidence level that all the y-ray 
flares occur at a low hard X-ray flux is > 7.7 normal standard 
deviations. 

We have also calculated the Discrete Cross-Correlation 
Function (DC F). The DCF corre lates two sets of unevenly 
sampled data (E delson et al. 198 "8*). We have correlated both 
BAT/Swift (sampled each day) and AMI-LA (variously sampled 
more than one time in a day or every few days) with AGILE data 
(sampled every few days, with a variable time bin size depend- 
ing by the y-ray state of the source). The calculation of the DCF 
uses a DCF binning of 8 days to take into account the biggest 
time bin size of AGILE light curve, with a scan time lag of 1 
day. The UL are converted into flux and related Icr error consid- 
ering half of the upper limit. The significance level of the peak 
was determined by simulating AGILE light curves with a boot- 
strap selection from the original points. 

Choosing the right DCF bin size is a trade off between the 
high accuracy of the correlation coefficient between sets of data 
of two time series and the resolution in the description of the 
cross-correlation curve. We have performed tests with different 
bin sizes of the AGILE light curve (splitting the same flux in 
smaller bins) and we have tried also different DCF bin sizes with 
the AGILE original light curve of Figure[3]and|5] In any case we 
have found that the peak of the DCF depends only weakly on the 
specific bin size. Due to the bin size of AGILE light curve, we 
use the DCF only for the determination of the cross-correlation 
and not for a measurement of a lag between hard-X and y-rays. 



3. Discussion 

The AGILE y-ray monitoring of Cyg X-3 contributes in a sig- 
nificant way to the determination of a repetitive pattern of y-ray 
emission and particle acceleration in this microquasar. Several 
fast (1- or 2-day) y-ray flaring episodes are detected in soft X- 
ray states. Furthermore, persistent and significant y-ray emis- 
sion is detected in AGILE data usually during a low flux in 
hard X-rays. The most intense y-ray emission above 100 MeV 
distinctively occurs in coincidence with minima of the hard X- 
ray emission. As we discuss below, on several occasions a ra- 
dio flare was detected near or after t he y-ray peak em ission. 
This trend, tha t was already apparent in dTavani et al. 2009b) and 
(lAbdo et al. 2009b) . is confirmed by a more extended monitor- 
ing over more than 1 year 

The y-ray flaring episodes are demonstrated to be even more 
relevant for understanding the source dynamics when compared 
to soft X-ray, hard X-ray, and radio data. Fig. [3] and |5] show the 
y-ray datapoints with multifrequency information. 

3. 1 . General characteristics 

Observationally, y-ray emission from Cyg X-3 starts to be de- 
tectable by AGILE when the Swift -BAT hard X-ray flux de- 
creases below a count rate flux of F - 0.02 cts. cm"^ s"', as 
shown by the data reported in Table |6] and Q (see Fig.[3]and|5j. 
Furthermore, y-ray flaring appears to be occurring during or im- 
mediately before sudden X-ray spectral transitions. The y-ray 
flare recorded on MJD 55324 is associated with a spectral tran- 
sition from hard to soft X-ray emission, although unlike the other 
detections, it is not coincident with a radio flare. 

The reported data show in general an anti-correlation be- 
tween the y and hard X-ray fluxes, and a correlation between 
the y and soft X-ray fluxeO 

Several y-ray flares (55001, 55019, 55034) are associated 
with 1-2 day delayed radio flares at the level of 1 Jy at 15 GHz 
as detected by AMI-LA. For other y-ray flares, the relation with 
the radio emission is not so obvious, either for lack of radio 
monitoring (see the episode at MJD 55025), or for the quasi- 
simultaneous y-ray and radio emission (as in MJD 55343, see 
Fig. 4). Ho wever, the majority o f Cyg X-3 y-ray flares reported 
here and in (iTavani et al. 2009bl) appear to follow a trend charac- 
terized by a substantial y-ray enhancement preceding radio flar- 
ing activity by a few days, which is within the formal DCF error 

We notice that occasionally during rapid variations of the X-ray 
flux the y-ray activity increases and becomes detectable at a -vT^ > 3 
even outside the soft X-ray state. See the case reported for MJD 55292- 
55300 and 55362-55370.. 



7 



A.Bulgarelli et al.: AGILE detection of Cygnus X-3 during the period mid-2009/mid-2010 



bar of 8 days. The y-ray flare on MJD 55019 is possibly associ- 
ated with a y-ray harder spectrum compared to other flares. 

An interesting case is provided by the event near MJD 
55025. A prominent radio flare reaching 4 Jy at 2.1 GHz 
(RATAN-600) and 1 Jy at 37 GHz (Metsahovi) occurred on MJD 
55030, followed by a secondary radio flare with inverted spec- 
trum near MJD 55035 detected by RATAN-600 as well as by 
AMI-LA and Metsahovi. AGILE data show that a 1-day y-ray 
flare was detected on MJD 55025, i.e., almost 5 days before 
the strong radio flare. Lack of radio monitoring preceding this 
flare does not allow us to determine the association of the MJD 
55025 y-ray flare with other radio flaring activity. In any case, 
during the ~4 days between the y-ray and radio flares the overall 
y-ray flux stayed at a low level. In particular, during the radio 
flare itself that lasted about 3 days, the y-ray flux was below 
F = 100 X lO"*^ ph.cm"^ s"' near 100 MeV and above. 

A sequence of y-ray emission immediately followed. On 
MJD 55034-55035 the y-ray flux shows an increase above F - 
180 X 10 ph.cm^^s"' for E>100 MeV followed by 4 days 
of high level flux and subsequent decay. One day later (MJD 
55035 and following days), the "secondary" radio flare with the 
inverted spectrum was detected reaching 2 Jy at 1 1.2 GHz. The 
excellent multi-frequency coverage by our group during this pe- 
riod (MJD 55030-55040) is enhanced by several TeV observa- 
tions of Cyg X-3 by the MAGIC group (A leksicet al. 2010 ) fol- 
lowed after an AGILE alert. It is then interesting to note that the 
MAGIC 95% confidence upper Umit of 4 x 10 '^ph.cm^^s ' 
above 250 GeV obtained from a set of observations including 
those of interest here provide, for the first time for Cyg X-3, si- 
multaneous broad-band spectral information from 100 MeV to 
TeV energies. 

The AMI-LA radio data of MJD 55343 (see Fig.© shows the 
correlation between radio and y-ray flares with the radio flares 
that anticipate the y-r ay flare. The same behaviour is reported in 
dWilliams et al. 201 lb . 

3.2. The y-ray spectrum 

Fig. |4] shows the y-ray spectrum obtained by summing all the 
major aforementioned episodes during the period 2009 Jun.-Jul. 
The complex relation between the y-ray emission of Cyg X-3 
near 100 MeV up to 1 GeV and the X-ray and radio spectrum 
will be addressed in detail elsewhere. We emphasize here that the 
AGILE results set interesting constraints on the hardness of the 
y-ray spectrum near 100 MeV. Our measured power-law spec- 
trum is consistent in being "flat", that is with a power-law index 
near 2 in the AGILE-GRID energy range. However, we cannot 
exclude the existence of spectral curvature of the vFy spectrum 
in the energy range below or near 100 MeV. We address this 
point as well as the analysis and i mplications of the m easured 
broad-band spectrum of Cyg X-3 in dPiano et al. in preph . 

3.3. AGILE and Fermi comparison 

We compare our results on Cyg X-3 with those of ferm/-LAT. 
There is considerable overlap of emission episodes detected by 
the two missions during the period 2009 Jun.-Jul. Fig. [3] (top 
panel) shows the AGILE y-ray lightcurves of Cyg X-3 with dif- 
ferent time bins an d the available Ferm;-LAT y-ray lightcurve 
(lAbdo et al. 2 009b). The two instruments have a quite different 
response and daily exposure at energies near 100 MeV, and the 
daily effective exposure is influenced by solar panel and other 
geometric constraints of the pointing strategy. Furthermore, for 



a source of relatively rapidly varying (within one day or shorter 
timescale) y-ray emission such as Cyg X-3, the two instru- 
ments can catch different states of emission. The 1-day flaring 
AGILE-GRID lightcurve generally agrees with the 4-day aver- 
aged ferwi-LAT lightcurve. Panel 2 of Fig. [3] shows the MJD 
55019 hard y-ray episode detected by AGILE (associated with 
the radio flare of MJD 55021 and ASM X-ray peak) which is not 
evident in the 4-day average of the corresponding Ferm/-LAT 
data shown in the top panel. 

Fig. ID (top panel) shows the time evolution of the off-axis 
angle of Cyg X-3 with respect to the instrument boresights of 
the AGILE-GRID (in red) and Fenni-L AT (in blue) for the 10- 
day interval covering the period 18-28 July, 2009 with respect 
to Cyg X-3 position for energies between 100 MeV and 400 
MeV. During this interval the position of Cyg X-3 was sta- 
ble near 30 degree off'-axis in the AGILE data. The Fermi-LAI 
pointing strategy is very different from the AGILE fixed point- 
ing mode strategy, and the ferm/-LAT off-axis angle changes 
continuously and rapidly for a typical 7-fold sampling of the 
source every day with a continuous sweep of the pointing di- 
rection. This difference in source pointing for Cyg X-3 in July 
2009 is clearly shown in the bottom panel of Fig. |7] where the 
cumulative effective area of AGILE-GRID is greater than the 
Fenni-LAiT (we have used the Pass 6 Version 3 Front photons 
Instrument Response Function). The cumulative exposure of the 
two instruments operating with the same pointing strategy with 
respect to Cyg X-3 position is still comparable (AGILE-GRID 
effective area is about half that of the Fermi-LAI). Fig.[8]shows 
the case where AGILE operates in spinning mode. 

Both the AGILE-GRID average flux above 100 MeV and 
spectral index are in good agreement with those reported in 
the Fermi-LAT First Catalog dAbdo e t al. 2010) but difficult to 
reconcile with the r esults published in ( Abdo e t al. 2009b) and 
(ICorbel et al. 2010ah who report an average flux above 100 
MeV of F ^ 50 x lO^^ph. cm"' s"' at the location of Cygnus X-3 
outside the active y-ray period, consistent with a stable emission. 
This relatively large flux would have produced a much stronger 
stable y-ray source in the AGILE-GRID data of the Cygnus re- 
gion. 

4. Conclusions 

Our present study of the y-ray emission from the microquasar 
Cyg X-3 adds substantially to the info rmation already gath ered 
from the first results of AGILE -GRID dTavani et al. 2009bl) and 
Fermi-LAI dAbdo et al. 2009bl) . We confirm the time variable 
nature of the extreme particle acceleration in the microquasar 
Cyg X-3 that manifests itself as y-ray emission above 100 MeV. 

The overall anticorrelation between the active y-ray states 
of Cyg X-3 and its hard X-ray e mission is evident i n our data 
that confirm the conclusions of dTavani et al. 2009"bl) . Both the 
AGILE pointing and spinning mode observations of the Cygnus 
region sampled Cyg X-3 in different spectral states, both "hard", 
and "soft". Enhanced y-ray emission is definitely produced only 
during the "soft" spectral states. These states are usually charac- 
terized by low values of the radio emission before major radio 
outbursts occur. 

We detected several y-ray flaring episodes from Cyg X-3 
with fluxes above F = 100 x 10"'^ ph. cm"^ s"' for photons en- 
ergies larger than 100 MeV. All these episodes occurred during 
a low flux in hard X-rays, and were simultaneous with either 
peak values of the X-ray intensity or low values of the radio 
flux that typically precede major plasmoid ejections and radio 
flares. Several episodes of enhanced y-ray emission with radio 
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flares were detected (MJD 54998, 55001, 55003, 55007, 55019, 
55034, 55343 and possibly 55025). In the majority of cases, (ex- 
cept MJD 55323 that is seems correlated with a fast transition be- 
tween X-ray states) the y-ray emission above 100 MeV precedes 
or is near the major radio flare during a rapid spectral change 
transition. 

The AGILE monitoring of Cyg X-3 in the spinning mode 
produced a consistent database with no significant data gaps 
since 2009 Nov. We have continuously sampled many X-ray 
state changes. Cyg X-3 stayed in the hard X-ray state most of 
the time since 2009 Nov. (see Fig.|5]l, implying a non-detectable 
y-ray emission from Cyg X-3 by AGILE. We finally stress that 
the AGILE monitoring without temporal gap is crucial for the 
comprehension of these transient phenomena. 

For the first time (around MJD 55000), we witnessed the 
transition from high to low flux in hard X-ray, with the conse- 
quent "ignition" of the y-ray emission process. Whether this is 
really an "ignition" , i.e., an enhanced and very efficient parti- 
cle acceleration process that start with the y-ray production, or 
a particle and photon transport change of the medium surround- 
ing the Cyg X-3 compact source, is a matter of debate. AGILE 
will continue to study Cyg X-3 in order to discriminate between 
theoretical models. 

5. Acknowledgements 

The AGILE Mission is funded by the Italian Space Agency 
(ASI) with scientific and programmatic participation by the 
Italian Institute of Astrophysics (INAF) and the Italian Institute 
of Nuclear Physics (INFN). Research partially supported 
through the ASI grants no. 1/089/06/2 and 1/042/10/0. We ac- 
knowledge financial contribution from the agreement ASI-INAF 
1/009/10/0. The Metsahovi team acknowledges the support from 
the Academy of Finland for our observing projects (numbers 
212656, 210338, and others). 



References 



Abdo, A.A., et al. (Fermi Collaboration) 2009, Science, 325, 840 

Abdo, A.A., et al. {Fermi Collaboration) 2009, Science, 326, 1512 

Abdo, A.A. et al. (Fermi Collaboration) 2010, ApJS 188, 405 

Aleksic J. et al. (MAGIC Collaboration) 2010, ApJ, 721, 843-855 

Argan, A., et al. 2004, Proc. lEE-NSS, 1, 371 

Barbiellini, G., et al. 2002, Nucl. Instr and Meth. A, 490, 146 

Benjamini, Y., Hochberg, Y. 1995, J. R. Stat. Soc. B, 57, 289 

Bulgarelli, A., et al. 2008, ASP Conference Series, 411, 362 

Bulgarelli, A., et al. 2010, Nucl. Instr and Meth. A, 614, 213 

BulgarelU, A., et al., 201 1 submitted. 

BulgarelU, A., et al. 2010, ATel n. 2609 

BulgarelU, A., et al. 2010, ATel n. 2645 

Camilo, F., et al. 2009, ApJ, 705,1 (C09) 

Cattaneo, P W., et al. 2011, Proc. of the RICAP 2009, 251 

Chen, A. W., et al. 201 1, A&A, 525 

Chen, A. W., et al. 201 lb, in prep. 

Chi X., Wolfendale A.W. 1991, J.Phys.G 17,987 

Corbel, S., et al. 2010, ATel n. 261 1 

Corbel, S. and Hays, E. 2010, ATel n. 2646 

Dame T.M. et al. 2001, ApJ, 547, 792 

Edelson, R. A., Krolik, J. H. 1988, ApJ, 333, 646 

Feroci, M., et al. 2007, Nucl. Instr and Meth. A, 581, 728 

Giacconi, R., et al. 1967, ApJ, 148, LI 19 

Grove, E., et al. 2008, ATel n. 1850 

Halpern, J.P, et al. 2008, ApJ, 688, L33 

Hjalmarsdotter, L., Zdziarski, A.A., Szostek, A., Hannikainen, D.C. 2009, 
MNRAS, 392, 251 

Kalberla, P M. W. , et al. 2005, A&A , 440:775-782 

Koljonen, K.I.I., Hannikainen, D.C, McCollough, M.L., Pooley, G., Trushkin, 
S.A. 2010, MNRAS, 406, 307 

Labanti, C, et al. 2006, Nuclear Physics B Proc. Suppl., 150, 34 

Giuliani, A., et al. 2004, Memorie delta Societa Astronomica Italiana 

Supplementi, 5, 135 

Giuliani, A., et al. 2006, Nucl. Instr and Meth. A, 568, 692 
Hopkins, A.M., et al. 2002, AJ, 123, 1086 

Mason, K.O., Cordova, FA., and White, N.E. 1986, ApJ, 309, 700 

Mattox, J. R., et al. 1996, ApJ, 461, 396 

Mattox, J. R., Hartman, R. C, Reimer, O. 2001, ApJ, 135, 155 

McCollough, M.L., et al. 1999, ApJ, 517, 951 

McLaughlin, M. A. , et al. 1996, ApJ, 473(2):763 

Miller C. J. et al. 2001, ApJ, 122, 349 

Mioduszewski A.J., Rupen, M.P, Hjellming, R.M., Pooley, G.G. 2001, ApJ, 
553, 766 

Molnar, L.A., Reid, M.J., and Grindlay, J.E. 1988, ApJ, 331, 494 

Parsignault, D.R., Gurski, H., Kellog, E.M., et al. 1972, Nature, 239, 123 

Perotti, F, et al. 2006, Nucl. Instr and Meth. A, 556, 228 

Piano, G., et al., in prep. 

Pittori, C, et al. 2009, A&A , 506, 1563. 

Prest, M., et al. 2003, Nucl. Instr and Meth. A, 501, 280. 

Sabatini, S., et al. 2010, ApJ, 712, L10-L15. 

Sanford, P W. and Hawkins, F H. 1972, Nature Phys Sci, 239, 135 

Szostek, A., Zdziarski, A.A., McCollough, M.L. 2008, MNRAS, 388, 1001 

Tavani, M., et al. 2009a, A&A , 520, 995 

Tavani, M., et al. 2009b, Nature, 462, 620 

van Kerkwijk, M.H., Charles, PA., Geballe, T.R., et al. 1992, Nature, 355, 703 
van Kerkwijk, M.H., Geballe, T.R., King, D.L., van der Klis, M., van Paradijs, 
J. 1996, A&A, 314, 521 

Vilhu, O. Hakala, P, Hannikainen, D., McCollough, M., & Koljonen, K. 2009, 
A&A, 501, 679 

Williams, P K. G., et al. 2011, ApJ, 733:L20 

Zdziarski, A.A. & Gierlinski, M. 2004, Prog. Theor Phys. Suppl., 155, 99 



9 



A.Bulgarelli et al.: AGILE detection of Cygnus X-3 during the period mid-2009/imd-2010 



Table 4. y-ray pulsars used to evaluate the instrument stability for E 100 MeV. 



AGL source name 


(l,b) positioning 






y-ray flux 




Counterpart name 




(degrees) 






(10 


ph. cm - s ') 






AGL J0247+6027 


(136.71, 0.70) ± 


0.37 = 


tO.l 


11 : 


t3 


4.4 


PSR J0248-h6021 


AGL J0303+7438 


(131.46, 13.98): 


tO.43 


±0.1 


7 ± 


2 


4.0 


PSR J0308-h7442 


AGL Jl 136-6053 


(293.90, 0.67) ± 


0.35 = 


tO.l 


12: 


t3 


4.5 


LAT PSR Jl 135-6055 


AGL J 1435-5932 


(315.68, 0.75) ± 


0.70 = 


tO.l 


13 : 


t4 


4.0 


LAT PSR J 1429-59 11 


AGL J1953+3254 


(68.76, 2.86) ± 0.18 ± 


0.1 


17 : 


t 2 


8.3 


PSR J1952+3252 


AGL J2224-I-6113 


(106.20, 3.20) ± 


0.43 = 


tO.l 


11 : 


t 2 


5.7 


PSR J2229+6114 



Table 5. 1-day y-ray flares consistent with the Cyg X-3 position reported in iTavani et al. 20699il for E>100 MeV. 



y-ray flaring date 


(l,b) positioninj 






y-ray 


flux 


FDR-ff 




(MJD) 


(degrees) 






(10-** 


ph. cm"^ s"') 






54572.00-54573.00 


(79.10, 0.60) ± 


0.60 H 


= 0.1 


260 + 


: 80 


0.001 


4.2 


54772.00-54773.00 


(79.30, 0.70) ± 


0.70 H 


= 0.1 


258 + 


: 82 


0.01 


4.0 


54811.00-54821.00 


(79.60, 0.30) ± 


0.60 H 


= 0.1 


210 + 


:73 


0.01 


3.8 



Table 6. Main 1-day and 2-days y-ray flares consistent with the Cyg X-3 position detected in 2009 Jun-Jul and from 2009 Dec. to 
2010 Junfor E>100 MeV. The last row reports the sum of all the flares with the repeated flare occurrence post trial signiflcance. 



y-ray flaring date 


(l,b) positioning 




y-ray 


flux 


FDR-ff 




pre-trial a 


post-trialo cr 


Swift -BAT flux 


(MJD) 


(degrees) 




(10-8 


ph. cm"^ s"') 






for a single 
detection 


for a single 
detection 


ictscm"^ s^^) 


54997.50-54998.50 


(80.50, 0.53) ± 0.90 


±0.1 


180 ± 


: 64 


0.005 


3.9 


3.75 


2.2 


0.010 


55001.00-55002.00 


(79.44, 0.93) + 0.65 


±0.1 


168 ± 


:67 


10-5 


3.5 


3.35 


1.57 


0.015 


55003.00-55004.00 


(80.35, 1.15) ±0.64 


±0.1 


157 ± 


56 


10-5 


3.8 


3.71 


2.11 


0.007 


55007.00-55008.00 


(79.30, 0.81) ±0.57 


±0.1 


176 ± 


64 


10-5 


3.7 


3.51 


1.83 


0.005 


55025.00-55026.00 


(79.75, 1.15) ±0.58 


±0.1 


167 ± 


70 


0.006 


3.3 


3.20 


1.29 


0.008 


55034.00-55035.00 


(80.12, 0.96) ± 1.00 


±0.1 


168 ± 


67 


0.005 


3.4 


3.27 


1.43 


0.013 


55324.00-55326.00* 


(79.12, 0.91) ±0.83 


±0.1 


170 ± 


70 




3.3 


3.20 


1.29 


0.028 to 0.011 


55343.00-55345.00t 


(79.89, 0.71) ±0.86 


±0.1 


290 ± 


: 103 




3.7 


3.51 


1.83 


0.020 to 0.018 


















repeated flare 




















occurrence 




















post-trial cr 




Sum of above detections 


(79.69, 0.72) ± 0.30 


±0.1 


160 ± 


40 




6.2 




7.2 





* ATel 2609. 
t ATel 2645 

G We calculated the post-trial significance for 150 trials, which is equivalent of about 1 year of total exposure time. 



Table 7. Main 1-day and 2-days y-ray flares consistent with the Cyg X-3 position detected in 2009 Jun-Jul and from 2009 Dec. to 
2010 Jun. for E > 400 MeV. The last row reports the repeated flare occurrence post-trial signiflcance. 



y-ray flaring date 


(l,b) positioning 


y-ray flux 




pre-trial a 


post-trial tr 


Swift -BAT flux 


(MJD) 


(degrees) 


(10-** ph. cm-2 s-') 




for a single 
detection 


for a single 
detection 


(ctscm^~ s^^) 


54997.50-54998.50 


(80.42, 0.41) ±0.60 ±0.1 


67 ±29 


3.9 


3.75 


2.2 


0.010 


55000.50-55001.50 


(78.88, 0.56) ±0.92 ±0.1 


45 ± 24 


3.0 


2.99 


0.89 


0.015 


55018.50-55019.50 


(79.75, 0.55) ±0.60 ±0.1 


65 ± 33 


3.7 


3.51 


1.83 


0.004 


55343.50-55345.50 


(80.07, 0.56) ±0.70 ±0.1 


94 ±43 


3.5 


3.35 


1.57 


0.020 to 0.018 












repeated flare 














occurrence 














post-trial cr 




Sum of above detections 


(79.8, 0.5) ±0.31 ±0.1 


58 ± 18 


5.1 




3.8 
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Fig. 3. AGILE/GRID, AGILE/SA, Fermi/LAT, Swift/BAT, XTE/ASM, AMI-LA, RATAN-600 and Metsdhovi Radio Observatory data of 
Cyg X-3 during the uninterrupted 2-months period 2009 Jun.-Jul. ( Panel 1 (top)) the AGILE-GRID j-ray light curve with a variable 
window time to put in evidence the Cyg X-3 flares of 1- day timescale, the other data and upper-limits are determined with 1-2- 
or 4-day time intervals; the Fermi light curve of 4-days timescale. E>100 MeV. (Panel 2) the AGILE-GRID j-ray light curve for 
E>400 MeV and the coverage of MAGIC observation ( gray area ). ( Panel 3 ) the hard X-ray light curve as monitored by BAT ( 15-50 
keV) and by Super-AGILE (20-60 keV) with a daily timescale bin. (Panel 4) the soft X-ray light curve as monitored by XTE-ASM 
(1.3-12.1 keV) for a 1-day integration time bin; (Panel 5) the ASM hardness ratio (5 - l2keV)/(3 — 5keV) data. (Panel 6) AMI-LA 
radio flux monitoring of Cyg X-3 at 15 GHz and the Metsdhovi Radio Observatory data at 37 GHz. (Panel 7) RATAN-600 radiq 
telescope data at difl^erent frequencies. The AGILE/GRID y-ray upper limits are at the 2 — cr level, the flux error bars are 1 — cr 
values. 
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Fig. 4. The AGILE-GRID y-ray spectrum of Cyg X-3 obtained by integrating all flaring episodes in 2009 Jun.-Jul. A single power- 
law model fitting gives a photon spectral index is 2.02 + 0.28. The open circle indicates the energy channel not used fijr fitting. 
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Fig. 5. Simultaneous AGILE/GRID, Swift /BAT and radio monitoring data of Cyg X-3 during the uninterrupted 18-months in the 
period 2009 Dec.-201 1 May. (Top panel:} the AGILE-GRID y-ray lightcurve for E>100 MeV. (Second panel:) Swift /BAT lightcurve 
( in counts per second in the energy range 15-50 keV). (Third panel:) AMI-LA radio flux monitoring of Cyg X-3 at 15 GHz. 



13 



A.Bulgarelli et al.: AGILE detection of Cygnus X-3 during the period mid-2009/mid-2010 





Fig. 6. The y-ray intensity map of the sum of all the Cyg X-3 flares reported in this paper, in Galactic coordinates for E>100 MeV. 
Pixel size— 0.1° with 5-pixel Gaussian smoothing. Cyan contours: AGILE-GRID 95% confidence level of LAT PSR J2032+4127; 
Green continuous contour: the Cyg X-3 flares 95% confidence level (statistical and systematic errors) for E>100 MeV. Green 
dashed contour: the Cyg X-3 flares 95% confidence level ( statistical and systematic errors) for E>400 MeV. Black crosses: Fermi- 
LAT (1-year catalog) sources (crosses shown for contours too small to be visible). The green contour has been calculated with the 
multi-source likelihood analysis method described in the text, using persistent sources reported in Tabled Blue cross is the position 
of Cyg X-3. 
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Fig. 7. Different pointing strategies of Cyg X-3 for AGILE-GRID ( in red) and Fermi-LAT ( in blue ). Upper panel: time evolution of 
the instrument off-axis angle with respect to the fixed Cyg X-3 position during the period 2009 Jul 18-28 when the AGILE operates 
in pointing mode. The curve in blue show the off-axis angle evolution for the all-sky scanning pointing strategy of Fermi-LAT. The 
curve in red show the off-axis angle for the AGILE-GRID fixed pointing strategy in 2009 Jul. Bottom panel: time evolution of the 
cumulative exposure between 100 MeV and 400 MeV for Fermi-LAT (blue curve, using public data and Pass 6 Version 3 Front 
photons Instrument Response Function) and AGILE-GRID ( red curve) assuming a starting time on 2009 Jul 18. 
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Fig. 8. Same pointing strategies of Cyg X-3 for AGILE-GRID (in red) and Fermi-LAT (in blue). Upper panel: time evolution of 
the instrument off-axis angle with respect to the fixed Cyg X-3 position during the period 2009 Dec. 3-13, when AGILE started to 
operate in spinning mode.The curves show the off-axis angle evolution for the all-sky scanning pointing strategy of AGILE GRID 
(in red) and Fermi-LAT (in blue). Bottom panel: time evolution of the cumulative exposure between 100 MeV and 400 MeV for 
Fermi-LAT (blue curve, using public data and Pass 6 Version 3 Front photons Instrument Response Function) and AGILE-GRID 
( red curve) assuming a starting time on 2009 Dec. 3. 
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